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ABSTRACT 

We have measured the near-infrared colors and the fluxes of individual pixels in 68 galaxies common 
to the Spitzer Infrared Nearby Galaxies Survey and the Large Galaxy Atlas Survey. Each galaxy was 
separated into regions of increasingly red near-infrared colors. In the absence of dust extinction 
and other non-stellar emission, stellar populations are shown to have relatively constant NIR colors, 
independent of age. In regions of high star formation, the average intensity of pixels in red-excess 
regions (at 1.25 /im, 3.6 /im, 4.5 /im, 5.6 /im, 8.0 /im and 24 /im) scales linearly with the intrinsic 
intensity of Ha emission, and thus with the star-formation rate within the pixel. This suggests that 
most NIR-excess regions are not red because their light is being depleted by absorption. Instead, they 
are red because additional infrared light is being contributed by a process linked to star-formation. 
This is surprising because the shorter wavelength bands in our study (1.25 /im-5.6 /xm) do not probe 
emission from cold (10-20 K) and warm (50-100 K) dust associated with star-formation in molecular 
clouds. However, emission from hot dust (700-1000 K) and/or Polycyclic Aromatic Hydrocarbon 
molecules can explain the additional emission seen at the shorter wavelengths in our study. The 
contribution from hot dust and/or PAH emission at 2 /im-5 /im and PAH emission at 5.6 jam and 
8.0 /im scales linearly with warm dust emission at 24 /im and the intrinsic Ha emission. Since both are 
tied to the star-formation rate, our analysis shows that the NIR excess continuum emission and PAH 
emission at ~ 1 — 8 /im can be added to spectral energy distribution models in a very straight-forward 
way, by simply adding an additional component to the models that scales linearly with star-formation 
rate. 



1. INTRODUCTION 

The signatures of star formation are expressed across a 
galaxy's spectral energy distribution (SED). Ultraviolet 
and optical light trace the recently formed stellar pop- 
ulations, infrared light maps out luminous dust grains 
that re-emit this same starlight, and radio emission 
traces magnetic fields produced by short-lived, massive 
stars. How well these tracers can be used to determine a 
galaxy's star-formation rate depends on the galaxy's lu- 
minosity, dust content, nuclear activity, and metallicity. 
Adding to the complexity, the timescale of each tracer 
depends on the emission mechanism, hydrogen is ion- 
ized by only the most massive stars on timescales of 10 T 
years, while dust emission can be heated by lower mass, 
longer-lived stars on the order of several 10 8 years. At 
a fundamental level these tracers are highly interdepen- 
dent, because energy must be conserved. For example, 
the energy in the UV radiation absorbed by dust must be 
balanced out by the re-emission of this energy at other 
wavelengths. This suggests that the best estimator for 
star-formation rate will be based on a combination of 
tracers which together capture the full energetics of the 
system. 



Rec e nt wo rk by Calzetti et al. (2007) and Kennicutt 



et al.| ( 2009 ) shows that a good estimator for the star 
lor mat ion rate emerges from the linear combination of a 
galaxy's near-UV or visible-wavelength emission (either 
UV continuum or line fluxes, attenuated by dust) and its 
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emission in the infrared (e.g. 8 /im, 24 /im, or total IRA 
8-1000 and also radio). The beauty of this approach is 
that the data itself can correct for the attenuation and 
emission of dust. 

At near-infrared (NIR) (0.9-5 /im) and mid-infrared 
(MIR) wavelengths (5-30 /im), emission due to heat- 
ing of very small grains (VSG) and polycyclic aro- 
matic hydrocarbons (PAH) also contributes to the spec- 
tral energy distribution. UV photons excite these 
small molecules which emit through various vibra- 
tional/stretching/rotational modes and thus this emis- 
sion should trace regions of star formation. PAH emis- 
sion around 8 /im shows spat ial correlation to sta r for- 
mation within H II regions ( |Calzetti et al. 2007), and 
can be used as a SFR indicator, particularly it used in 



combinatio n with an optical SFR indicator (Kennicutt 
et al.||2009| ). Observations of nearby star forming galax- 



ies suggest that the spectrum of PAH line emission is 
expressed through constant line ratios. For example 



|et al. ([2003]) found constant ratios of stellar subtracts 
tiux densities at 4, 6.2, 7.7 and 11.3. PAH and/or dust 
emission at shorter wavelengths has received less atten- 
tion, partly because of the stronger contribution from 
stellar continuum emission at these wavelengths. 

In addition to PAH line emission, NIR continuum emis- 
sion in excess of NIR stellar emission can be significant at 
wavelengths of 2-5 /im. High redshift star forming galax- 
ies, where SFRs are significantly higher than n earby, 



emit significant NIR excesses at this w avelength (Men- 



tuch et al. 2009 Magnelli et al. 2008) and its emi ssion 
corre lates with optical SFR indicators (iMentuch et al. 



2009). Nearby galaxies, studied with the Infrared Space 
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Observatory al so show NIR excess continuum emission 
( Lu et al.|2QQ3 ) which scales linearly with the PAH spec- 
trum and shows trends of increasing with increasing star 
formation indicators. This NIR dust emission component 
is seen not only in extragalactic emission, but has been 
found locally in galactic c irrus emission (iFlagey et al.| 
[2006] ), reflection n ebulae flSellgren et al f|1983[ 
planetary nebulae dde Ruyter et al. 2006 ) ^ana massive 
star-forming reg i ons| Maercker &; Bu rton| " 



cker et al. 



two cases, 



2006 
tEeN 



|Longmore et al.| |2007 ) . 
IK dust emission is lin] 



2005 Maer 



n the latter 
linked to circum- 



stellar disk emission (|Wood et aL||2008| |Touhami et al. 
20 1 51 jAcke fc van den Ancker|2004| where high luminos- 



ity OB stars heat dust to its sublimation temperature 
(for example a 10 5 L star will heat dust to 1000 K at r 
= 7.3 AU, while a 10 2 L star will heat dust to 1000 K at 
r = 0.23 AU and a L star will heat dust up to 1000 K 
out to a radius of 0.023 AU.) 

The source of extragalactic NIR dust emission, and its 
contribution to a galaxy's integrated light remains un- 
certain. The excess seen in higher star forming galaxies 
at high redshift is possibly a combination of both NIR 
dust emission related to star formation regions (which 
includes reflection nebulae and circumstellar dust emis- 
sion) and N IR emission due to cirrus emission. In Men-| 
tuch et al.| (|2009), we show that there is a trend ol in- 



creasing NIK continuum emission with star formation 
and show that the magnitude of NIR excess emission can 
be well matched if the majority of the excess NIR emis- 
sion originates in hot (700 K-1000 K) thermal emission, 
possibly from dust in circumstellar disks heated to their 
sublimation temperature. Other sources of NIR emis- 
sion such as dust heated by active galactic nuclei, plane- 
tary nebula and reflection nebula can also be important, 
although their contribution to the integrated light of a 



galaxy is expected to be less significant (Mentuch et al. 



2009]). 

7v central purpose of the present paper is to test this 
hypothesis, using the images of nearby galaxies in which 
we are able to resolve the components of the galaxy 
at greater resolution, giving us the opportunity to de- 
termine if the NIR excess emission is associated with 
emission from young stellar populations, galactic cirrus 
or active galactic nuclei. Additional narrow-band Ra 
imaging, which traces the young, massive stellar popula- 
tions, complements the analysis. We perform multiwave- 
length pixel-by-pixel analysis on a sample of 68 nearby 
galaxies fro m the Spitzer Infrared Nearby Galaxies Sur- 
vey (SINGS; |Kennicutt et al.|2 003) with complementary 
near- infrared' data from the Large Galaxy Atlas (LGA; 
jJarrett et al.|2003| ). These data resolve galaxies on phys- 
ical scales between 8 pc (for the closet galaxy NGC 6822 
at a distance of 0.6 Mpc) and ~ 400 pc (for the most 
distant galaxies at ~30Mpc). 

Near-infrared colors of evolved stellar populations are 
not highly sensitive to a population's age, although con- 
tributions from post-main-sequence stars can contribute 
at ages of ~l-3Gyr. The NIR traces the Rayleigh- Jeans 
tail of stellar emission and as a result the colors remain 
nearly constant as stars of progressively lower masses 
leave the main sequence. However, observations have 
shown that in addition to emission in the NIR, extinc- 
tion due to dust leads to reddened colors as is shown 



in Indebetouw et ah fl2005|) in our ow n galaxy and in 
nearby starbursts ( jCalzetti et al.||1994 ). Integrated col- 
ors consist of light from all the stellar populations and 
dust in a galaxy, making it difficult to discern whether 
red NIR colors are associated with non-stellar emission 
or are merely due to dust extinction. Our study provides 
synergy between investigations of integrated emission for 
whole galaxies, and studies of resolved stellar populations 
by grouping pixels with common near-infrared colors to- 
gether so that the common properties of these regions 
are revealed. 

We motivate the need for pixel-by-pixel color studies 
in the following section. The rest of the paper is out- 
lined as follows. In §[3j we briefly summarize our sample 
selection from the SINGS and LGA surveys. Our pixel- 
by-pixel color selection and image processing is described 
in §[4j Results from the average NIR colors and intensi- 
ties at Ha, 1.25 /im, 3.6 /im, 4.5 /im, 5.6 /im, 8.0 jam and 
24 /im measured in each color selected region are pre- 
sented in §[5] Implications for studying integrated galax- 
ies and applying our results to SED models are discussed 
in §[6} 

2. PIXEL-BY-PIXEL NIR COLORS OF RESOLVED 
GALAXIES 

An analysis of galaxy colors at the level of individual 
pixels is substantially more complicated than the usual 
approach based on integrated photometry using large 
apertures. It is important to establish from the outset 
whether such a complicated approach is really required 
in order to understand the underlying physics. Before 
embarking on a description of our methodology, we will 
first attempt to justify the underlying need for it. 

At near-ultraviolet and visible wavelengths a galaxy's 
spectral energy distribution is the sum of the light from 
the stellar populations that make up the galaxy, modu- 
lated by dust absorption, and augmented by additional 
sources of emission (such as from nebular lines). The 
different ages of the stellar populations impact the SED 
and the resulting colors, as young high-mass stars peak 
in their emission at shorter wavelengths and have blue 
visible colors, while older low-mass stars peak in their 
emission at red wavelengths and have redder visible col- 
ors. Slightly longer near-infrared wavelengths (around 
~l-2 /im) are ideal for probing the bulk of a galaxy's 
stellar mass, because starlight still dominates the overall 
SED and most of the light at these wavelengths comes 
from late evolutionary stages of low-mass stars which 
dominate the stellar mass budget. The colors of stel- 
lar populations in the near-infrared are relatively con- 
stant for most stellar spectral types and most evolution- 
ary phases , with few notab le exceptions like the TP-AGB 
stars (e.g. |Maraston| [2005 ) . Slightly longer wavelengths, 
say between 2 /im and 5 /im, are special. This wavelength 
range corresponds to a minimum in a galaxy's SED, 
because the Rayleigh- Jeans tail of starlight is rapidly 
declining but thermal emission from warm dust is not 
yet important. Additional components of a galaxy that 
happen to emit at these wavelengths, such as hot dust, 
may be detectable at these wavelengths even if they do 
not contribute a large fraction of a galaxy's bolomet- 
ric output, simply because of contrast against the weak 
underlying galactic SED. At wavelengths longer than 
5 /im starlight starts to play a negligible role, and the 
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physics of dust emission takes over. For example in 
the Spitzer 5.6 fim and 8.0 fim bands, molecular emission 
lines from Polycyclic Aromatic Hydrocarbons (PAHs) are 
important, while at 24 /am most of the light comes from 
the thermal emission of warm dust. 

Some of these ideas are illustrated in Figure [I] The 
left-hand panel of the figure shows a ratio of optical and 



near- infrared colors. Integrated photometry from |Dale 
et al. (2007) for the SINGS sample (black asterisks) is 



plotted along with rest-frame colors from a high-redshift 
galaxy sample (data taken f rom the Gemini Deep Deep 
Survey; [Abraham et al.|2004[ ). For the high-redshift data, 
the colors were compu ted from the best fit SEDs from 
Mentuch et a l. (2009) and have been corrected for full 
cosmological dimming. Their symbols have been keyed 
to differing galax y types as indicated by the spectral in- 
dicators given in [Abraham et ah] ( 2004 ) . Star forming 
galaxies show spectral indicators of young stellar popu- 
lations, while evolved galaxies are mainly comprised of 
older stellar populations. Intermediate galaxies show sig- 
natures of both evolved and young stellar populations. 
The axes can be converted to AB magnitude colors by 
multiplying by 2.5. 

The spread in the I1.25/I0.43 ratio (or, in magnitude 
space, B — J color) is over 1 dex and is much larger than 
the spread in the NIR intensity ratios (J — H color). 
Star forming galaxies tend to have bluer B — J color and 
evolved galaxies have redder B — J color, but the inter- 
mediate galaxies are spread throughout (dust can play 
a major role in this spread). In the right-hand panel, 
near-infrared colors are plotted for both the SINGS and 
GDDS rest-frame sample. The rest-frame colors of the 
high redshift data are interpolated from SED models fit 
to the observed photometry. The purple crosses are from 
SED models that only incorporate emission from star 
light and attenuation by dust. The NIR colors of these 
models colors occupy a compact region i n the NIR color- 
color plot. However, in Mentuch et al. (2009), we show 
that these models are not a good description of the near- 
infrared observed photometry and find that a SED stellar 
population augmented with a NIR emission model better 
represents the observations. Rest-frame colors derived 
from these SEDs are shown in the right panel (keyed 
to galaxy spectral types). These rest-frame colors of 
the high-redshift GDDS sample are consistent with the 
SINGS photometry for galaxies composed primarily of 
evolved populations, but for the star forming and inter- 
mediate GDDS galaxies, which have higher SFRs com- 
pared to the low redshift SINGS dataset, show a large 
spread toward increasing I3.6/I1.25 ratios (or bluer J — L 
color) and I4.5/I3.6 ratios. 

On the whole, Figure [T] makes three main points which 
motivate the analysis in this paper. (1) The purely stellar 
component of the SED (derived from fitting of SPS mod- 
els, purple crosses in the right panel) shows a small (~0.2 
dex) variation in the NIR which is much smaller than at 
visible wavelengths, where the variation is ~1 dex. (2) 
In spite of the previous point, real galaxies show much 
greater variation (~ 1 dex) in near-infrared color than 
can be encompassed by stellar population models. (3) 
As a result of the previous two points, the large variation 
in the NIR galaxy colors must be telling us quite a lot 
about dust or some other source of non- stellar emission. 
Either the starlight is being reddened by dust absorp- 



tion, or dust emission must be contributing to the colors 
at surprisingly short wavelengths. 

Resolved galaxy colors can help distinguish between 
the latter two possibilities, because the properties of dust 
are known to vary from place-to-place in galaxies. Fig- 
ure fl] (right panel) shows that galaxy colors tend to scat- 
ter along the extinction vectors found in starburst re- 
gions, but not along the vectors appropriate for regions 
outside of starburst regions. The extinction curve for 
lines of sight thro ugh our own galaxy (ob tained by the 
GLIMPSE survey, |Indebetouw et al7||2QQ5| ) is shown in 
green. Using the prescriptions presented in |Pei| ([1992), 
we see that this extinction curve is consistent with the 
dust extinction law for the Milky Way (shown in red), the 
Small Magellanic Cloud (shown in orange) and the Large 
Magellanic Cloud (shown cyan). The extinction mea- 



sured in starburst galaxies from Calzetti et al. (1994) is 
plotted in blue, and th e effective dust extinction law from 



Chariot fc Fall (2000) is shown in black. In this figure 
we show all extinction vectors originating from the colors 
(logio I3.6/I1.25, logio I4.5/I3.6 = -0-30, -0.17) which we 
found to be the average colors of stellar emission in the 
NIR from the PEGASE.2 stellar SED^l The trend for 
NIR colors along this vector, predicted to be due to dust 
extinction in starburst regions, is supported by the high 
redshift data. The star forming galaxies tend to have the 
reddest NIR colors, but we can not say for certain that 
dust extinction is the main culprit. Most of the GDDS 
galaxies have total extinction values A v < 1.5 mag. The 
extinction vectors here suggest the extinction is in excess 
of A v > 2 mag. What do the resolved colors of nearby 
galaxies tell us about dust extinction in the NIR? 



3. SAMPLE SELECTION 

Table [1] summarizes our sample of 68 spatially resolved 
nearby galaxies, with observations at 1.25 /im (J-band), 
3.6 /im, 4.5 /im, 5.6 /im, 8.0 /im and 24 /im. All data come 
from publicly available sources, namel y the Spitzer In- 
frared Nearby Galaxies Survey (SINGS; |Kennicutt et al. 
|2003[ ) and the Two Micro n All Sky Survey (2M ASS) 
Large Galaxy Atlas survey ( Jarrett et al.|2003| hereafter 
referred to as the 2MASS LGA). We used SINGS data 
release 2070410_enhanced_vl4. Our data is comprised of 
all galaxies that were observed by both surveys that meet 
our minimum signal-to- noise criterion (described below). 
As a result we excluded SINGS galaxies DDO 053, Holm- 
berg I, Holmberg IX, M81 Dwarf A and M81 Dwarf B, 
which do not have NIR coverage in the 2MASS LGA. In 
addition, two SINGS galaxies, DDO 154 and DDO 165, 
did not contain any pixels which satisfied our signal-to- 
noise criterion, so only two objects were omitted for our 
analysis because of low signal-to- noise. 

In addition to the infrared data contained in the 
SINGS and 2MASS LGA data releases, 56 galaxies had 
publicly-available Ka images, obtained from telescopes 
at Kitt Peak National Observatory (KPNO) or Cerro 



3 These colors are likely redder than that from pure stellar emis- 
sion. Attenuation due to dust was included in these models. We 
later show (in §[5j that these colors are consistent with the median 
pixel color of the majority of pixels which compose the SINGS 
galaxy sample. 
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a Shows spectral signatures of an active galactic nuclei ( Dale et al.||2006| 
b Indicates Ha observations 
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Fig. 1. — Color-color plots in the optical/NIR (left) and the NIR (right) for the integrated SINGS photometry from pale et al.N2007t 
and r est-frame colors for the GDDS high-redshift galaxy sample. Rest-frame colors are computed using best fit SEDs from Mentuch et al. 
(2009). The spread in I1.25/I0.43 color (left panel) results from variations in both age and dust attenuation and is well described by SPS 
models, while the spread in I3.6/I1.25 color (right panel) is largest for high SFR galaxies, but is not well described by SPS models (purple 
points) although SPS models augmented with NIR emission (colored points, keyed to legend according to galaxy spectral type) provide a 
better model. In the right panel, extinction vectors are plotted from the color locus of (logio I3.6/I1.25, logio I4.5/I3.6 — -0.30, -0.17) for 
various dust extinction laws. The colored boxes indicate regions of increasing NIR color, defined in §^]and used as color selection in the 
quantitative analysis presented in this analysis. 



6 



TABLE 2 

Color key for Near-infrared color 
criteria 



Color Key 


a log 10 I3.6/I1.25 b lc 


)gio I4.5/I3.6 


Red 


-0.6 to -0.4 


-0.19 


Orange 


-0.4 to -0.2 


-0.17 


Yellow 


-0.2 to 0.0 


-0.14 


Green 


0.0 to 0.2 


-0.12 


Blue 


>0.2 


-0.09 


Purple 


Manually selected 





a Upper and lower limits. 

b Value in the center of the parallelogram shown 
in Figure fl] 

Tololo Inter- American Observatory (CTIO). We used 
the continuum subtracted Ha images contained within 
the SINGS data release, applying a small correction for 
a residual background pedestal where needed. 

Because nuclear activity may contribute significantly 
to NIR and MIR flux, we took careful consideration of 
which systems hosted active galactic nuclei. Galaxies 
were differen tiated based on n uclear activity using the 

20061), who used optical line 
SINGS nuclei are either 



results from Dale et al 
diagnostics to identity wmcl 
starburst dominated or require additional powering from 
an AGN, using the Liner /Starburst boundary of N II 
A6583/Ha ~ 0.6. This identification is included in Ta- 
bled] 

4. METHODS 

As we have noted in §[2] and illustrated in Figure [I] the 
scatter in the NIR colors of high-redshift star forming 
galaxies is greater than that seen nearby (where SFRs 
are lower). How does this small scatter manifest itself 
on smaller spatial scales? Our methodology is to group 
the pixels composing the local galaxy sample into re- 
gions based on different color-cuts motivated by both the 
spread in the data and the extinction law from Chariot 
fc Fall ( 2000 ) as indicated by the colored parallelograms 
shown m Figure [I] The upper and lower limits of the 
boxes are defined as 0.1 dex ab ove an d below the extinc- 

( [2000] ). 



tion vector of jCharlot fc Fall 
separated in steps ol increasing logio I3.6/I1.25 = 
(or 0.5 mags in L — J color) from logio I3.6/I1.25 



The regions are 
0.2 dex 
= -0.6. 



We summarize the properties of the NIR color-selected 
regions in Table [2] 

4.1. Image preparation 

Prior to analyzing the pixel-by-pixel colors of our 
galaxy sample a number of steps were taken to regis- 
ter images to a common pixel scale and common point- 
spread- function (PSF), clean out contaminants and en- 
sure that each pixel has a signal-to- noise (S/N) ratio of 
at least 15. Each step is outlined below in more detail. 

1. All images were re-sampled to the IRAC plate 
scale of 1.22 "/pi xel using th e remap tool in the 
WCStools library (Mink 2002]). 



2. The counts in each pixel were converted to intensity 
units of MJy/sr. 

3. Images were registered to sub-pixel accuracy us- 
ing unsaturated bright stars and 500 by 500 pixel 



(^10.2 1 by ~10.2') subsets were cut from the reg- 
istered master frames. 

4. Saturated regions and large galaxy neighbours on 
each image were manually masked. 

5. All images were convolved with a gaussian ker- 
nel to match the seeing-limited ~6" beam of the 
24 /im MIPS image point-spread function. We as- 
sumed a median FWHM of 1.0", 2.8", 1.7", 1.7", 
1.9", 2.0" and 6.0" for the 1.25, 3.6, 4.5, 5.6, 8.0 
and 24 /im images, respectively. (Pixel- by-pixel im- 
age analysis requires that each pixel value be a good 
representation of the true intensity in the area of 
galaxy contained by the pixel. By convolving our 
images, we ensured that the distribution of inten- 
sities were matched as they were spread across the 
pixels by varying point-spread functions.) 

6. Bright foreground stars and small background 
galaxies were automatically removed. Foreground 
stars were removed in 5 pixel (6.1") radius aper- 
tures using the the 2MASS point source cata- 
log. Unfortunately many compact H II regions are 
falsely identified as point sources in the 2MASS 
catalog, so we attempted to identify these sources 
beforehand through their m id-infrared colors as in 



Muhoz-Mateos et al. 



_ ([2009]). We also did not re- 
move any point sources with 24 /im intensity > 
0.7 MJy. These two criteria result in ~30% cata- 
log point source removal, with some contaminants 
remaining, but we chose this approach to preserve 
compact H II regions and compact nuclear sources 
such as AGN. In addition, faint background galax- 
ies have the potential to show up as spurious ob- 
jects in our analysis. We have inspected each im- 
age carefully to note cases where contamination re- 
mains important after the steps just described, and 
will explicitly note these cases below. 

7. We used a m edian adaptive smoothing algorithm 
dZibetti|[2009l ) on the 1.25, 3.6, 4.5 /im images so 
that "SJly > 15 is obtained in each pixel. This final 
step is non-standard and will be described below. 

4.2. Adaptive smoothing 

Much of our analysis is based on an investigation of the 
flux ratios of individual pixels. This required the con- 
struction of data frames obtained by dividing 1.25 /im, 
3.6 /im, and 4.5 jam images against each other. To mini- 
mize the noise in the resul ting divided images we used the 
Adapt smooth software ( Zibetti|[2009[) to median smooth 
the images so that a reasonably nigh signal-to-noise 
(S/N) threshold is first met. This trades resolution at 
low signal-to-noise levels for robustness in the final flux 
ratios (see the discussion in Zibetti et al. 2009 for de- 
tails). 

Our goal with Adapt smooth was to obtain input data 
frames where S/N > 15 in each pixel. Adapt smooth 
works by taking the median intensity over a smoothing 
aperture with a radius defined by the size needed to meet 
the S/N > 15 criterion. If the smoothing radius required 
is greater than 12 pixels, the algorithm is declared to have 
failed for that pixel, and the pixel is excluded from sub- 
sequent analysis. All results reported in this analysis in 
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which colors were computed in individual pixels use the 
smoothed then convolved images. For analysis involving 
measurements of the average intensity over many pixels 
adaptsmoothed images are used to perform pixel selec- 
tion (which would be impossible otherwise) , whereas un- 
smoothed images are used to derive pixel intensities to 
be averaged. In fact, averaging over many pixels makes 
the adaptsmoothing step superfluous. Moreover, this ap- 
proach allows us to retain the most accurate flux infor- 
mation in each pixel and avoid possible biases (alth ough 
these biases are modest, as discussed in |Zibetti||2QQ9| ). 

5. RESULTS 

The near-infrared colors in the individual pixels of the 
68 galaxies in our sample were measured by performing 
image division with the median smoothed (and signal- 
to- noise optimized) then convolved images (as described 
in §[4|. The resulting pixel colors in the NIR color space 
are shown for 6 example galaxies from our sample in Fig- 
ure [2] Results for all galaxies can be found on the online 
edition of this article. We find 10 galaxies to have atyp- 
i cal p ixel distributions and these are discussed later in 



5.2 



In the leftmost column, 10.2' x 10.2' images of the 
galaxies are shown at IRAC 3.6 /im. Regions with satu- 
rated objects, foreground stars or neighbouring galaxies 
are removed. In the column to the right of the image col- 
umn, we indicate the spatial location of the NIR color- 
selected pixels. The colors are keyed according to the 
NIR color-selection, shown in Figure [I] and defined in 
Tabled 

The three rightmost columns show binned color-color 
diagrams with the intensity of cells keyed to various 
quantities. The intensity in the first of these three is 
keyed to the total number of pixels in that color bin. 
The scatter in pixel NIR colors is larger than for the 
integrated colors of the SINGS galaxies shown in Fig- 
ure [T] with some pixels containing intensity ratios at 
I3.6/T1.25 > 1, larger than any of the integrated colors. 
As with the integrated colors, the majority of pixels in 
the galaxies are found in the orange region, that is 0.1 dex 
(or 0.25 mags) above and below our typical color for the 
sample of logio I3.6A1.25j logio I4.5/I3.6 = -0-30, -0.17, or 
(1.25]-[3.6]ab = -0.75 and [3.6]-[4.5]ab = -0.425). For the 
spatially resolved galaxy sample, we find the majority of 
pixels are found at logio I3.6/I1.25 = —0.30 ± 0.07 and 
logio I4.5/I3.6 = —0.19 ± 0.02, consistent with the colors 
we visually selected from the integrated colors to repre- 
sent the colors of evolved stellar populations. A smaller 
subset of pixels tend to scatter along the extinc tion vec- 
tor defined by A -0 7 from empirical modeling by Chariot 
fc Fa ll (2000) , shown to be a goo d match to star burst 
galaxies from Calzetti et al.| ( |1994[ ). 

Figure [2] allows us to map the color-selected pixels into 
the spatial domain of the galaxy. The majority of pixels 
with NIR colors (keyed as orange regions) described by 
stellar emission are found throughout the galaxy. Re- 
gions with redder colors (identified as yellow then green 
then blue) follow the spiral structure and regions of star 
formation in the galaxy. A degeneracy between extinc- 
tion and emission occurs at the brightest knots of star 
formation, where the NIR colors are their reddest. If 
these regions correspond to regions of high dust extinc- 
tion, then the red NIR colors may be the result of dust 
extinction in which the light at 3.6 /im is attenuated less 



than the 1.25 /im emission. On the other hand, addi- 
tional emission can also explain the excess continuum 
emission at 2 — 5 /im. 

Fortunately, for 56 of the galaxies, continuum sub- 
tracted Ha images can be combined with 24 /im images to 
create two composite images that represent the intrinsic 
Ha emission (corrected for attenuation due to dust ex- 
tinction) and the dust extinction at Ha in each pixel. 
The intrinsic Ha line emission luminosity, and conse- 
quently the star formation rate, of normal star forming 
galaxies and H II regions can be estimated empirically as 
a linear superposition of the observed Ha a nd 24 /im flux 
([Kennicutt et al.|2007l [Calzetti et al.|2007[ ). Recent work 
([Kennicutt et al. 2009) shows that lor a sample of galax- 
ies, including the SINGS galaxies, the dust corrected, in- 
trinsic Hai nt luminosity is related to the observed Ho; 5 S 
and 24 jam luminosity by: 



L(Ha int ) = L(Ha obs ) + 0.02 x L(24/xm) (1) 



with the galaxies showing a rms dispersion of 0.119 dex 
in the fit. From this equation, one can immediately get 
Ina,int shown in Equation 2. Also, from L(H a ^ nt ) one 
can obtain t he SFR by mean s of a simple multiplicative 



factor (e.g. Kennicutt 19981). Hence, lua,int: is in fact 



a direct estimate ol the local star formation rate density 
(or surface brightness). 



lHa,obs + 0.02 X 



24.0 
0.656 



x I 



24 /im 



(2) 



Similarly, the average extinction in each pixel can be de- 
rived using the intrinsic Ha intensity, i"Ha,mt> defined 
above and the observed Ha intensity, In a ,obs • The ex- 
tinction, A(Ha) in mag, is given as: 



A{Ha) = 2.5 log 



10 



0.02 x (24.0/0.656)-?^ 

J-Hcx^obs _ 



(3) 



The two panels on the right of Figure [2] show the me- 
dian intrinsic Ha intensity, Ia a ,int^ an d the median dust 
extinction, A(Ha), in the pixels corresponding to the 
given NIR colors. By spatially resolving the galaxy col- 
ors, we can see a trend emerge in the NIR color-color 
plane in the pixel plots shown in Figure [2] The pix- 
els which scatter along the extinction line to the redder 
color-selected regions also tend to have the highest in- 
trinsic Ha intensities. Regions with higher intensity star 
formation show redder NIR colors. The dust extinction, 
in the far right panel, shows less variation, and tends to 
be A(Ha) =0 — 2 mag. There is no strong trend of 
increasingly red NIR colors with increasing dust extinc- 
tion. We now attempt to quantify this by measuring the 
average surface brightness in increasingly red NIR color 
regions. 
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Fig. 2. — Images at IRAC 3.6 /xm and the spatial locations of the NIR color selected pixels (defined in Table [2| are shown for 
examples of normal star-forming galaxies in our sample: NGC 0628, NGC 2976, NGC 3184, NGC 3627, NGC 519¥and NGC 6946. 
The three rightmost panels plot the pixel NIR colors with the intensity levels defined as the number of pixels at each color bin 
in the first of these three panels. NIR color selection criteria is shown in colored boxes. The last two panels show the median 
intrinsic Ka intensity (and consequently star formation rate density) and the median A (H a) in mag for all the pixels in that 
NIR color bin. Figures for all galaxies can be found in the online edition. 



9 



..Mr 



log 10 T H (MJy/sr) 



log10 T H (MJy/sr) 



A(Ha) 



A(Ha) 



Fig. 3. — The average I3.6/I1.25 (left column) and I4.5/I3.6 (right 
column) colors in each of the 5 regions denned in Figuref2]as a func- 
tion of the average dust-corrected Ha intensity (top panels) and 
dust extinction at Ha (bottom panels). The colors are keyed to the 
NIR color selected regions denned in§[4] and described in Table [2] 
The purple dots are regions with anomalous NIR color sequences 
as discussed in § |5.2| By design, the colors should change from red 
to blue with increasing I3.6/I1.25 and I4.5/I3.6 ratios. The dust 
extinction is scattered with values ranging from to 6 mags of 
extinction in all color-selected regions. There is a stronger corre- 
lation between redder NIR colors with the average dust-corrected 
Ha intensity than average dust extinction. 



5.1. Average surface brightness measurements in NIR 
color- selected regions 

We measured the average NIR colors and intensity 
(i.e. surface brightness) for each galaxy in each of the 
NIR color selected regions marking 0.1 dex in color space 
above and b elow t he extinction curve of A -0,7 from lChaM 
lot fc Fall| fl2000| ) and defined in Table [2] The regions 
are mapped in the spatial domain and in the NIR color 
plane in the second and third columns of Figure [2] Mea- 
surements are made at all wavelengths in this analysis 
(0.6563, 1.25, 3.6, 4.5, 5.6, 8.0 and 24 pm) as well as to 
the intrinsic Ha image defined by Equation [2] and the 
A(Ha) dust extinction image defined by Equation [3] 

We begin by looking at the NIR color trends. A trend 
of increasing I3.6/I1.25 and/or I4.5/I3.6 colors with in- 
creasing dust extinction would suggest that red NIR col- 
ors are due to dust, or alternatively, a trend of increas- 
ingly red NIR colors with the average dust-corrected Ha 
surface brightness would suggest that the NIR colors are 
instead due to dust emission. Because of the likely co- 
incidence of both high dust extinction and strong Ha 
emission in regions of star formation, it is not easy to dis- 
sociate these parameters. Our attempt to do so is shown 
in Figure [3J Here, the NIR colors in the color-selected 
regions are compared to the average dust-corrected Ha 
emission and dust extinction. The top row shows the 
NIR colors (I3.6/I1.25 an d I4.5/I3.6 from left to right) both 
increase with increasing Ha emission, as they do with 
dust extinction in the bottom rows. The data points 
are keyed to the color selection regions. As expected 



(by design) the regions with higher NIR ratios are from 
the redder color-selected regions (blue, green and yellow 
points). 

In all cases, there is a trend of increasing NIR col- 
ors with both Ha emission and dust extinction. We 
find the correlation is stronger between the NIR colors 
and the Ha emission than the dust extinction. Quan- 
tified by Pearson's correlation coefficient, the NIR in- 
tensity ratios, I3.6/I1.25 an d I4.5/I3.6? correlate to the 
intrinsic Ha surface brightness by p = 0.38 and 0.36, 
respectively. We compare how the NIR intensity ra- 
tios correlate with the dust attenuation, quantified as 
lHct int /lHct obs (to compare linear quantities), and find 
that the correlation is weaker, but still significant at 
p = 0.11 and 0.16 for I3.6/I1.25 an d I4.5/I3.6 ratios re- 
spectively. The stronger trend with Ha emission suggests 
the NIR colors are due to dust emission rather than ex- 
tinction. The weaker trend for the reddest color- selected 
regions to have a higher median extinction value is likely 
a consequence of the higher star formation and associated 
dust content in these regions, resulting in higher dust ex- 
tinction. 

The intensity in each region is compared for each ob- 
served waveband to the intrinsic Ha intensity (a direct 
measure of the average star formation rate) in Figure pi 
There is clear correlation between the intensity at 5.6, 8aJ 
and 24 pm and the intrinsic Ha intensity. These bands 
are dominated by PAH emission. Also evident is a corre- 
lation between emission at 1.25, 3.6 and 4.5 pm and the 
intrinsic Ha emission, although here the trends are more 
nuanced. In regions of high intrinsic Ha intensities (and 
high star formation rates), the emission from the young 
star-forming environment dominates over the emission of 
the longer lived population. But at lower Ha intensities 
(<10 MJy/sr), the emission at 3.6, 4.5 pm and particu- 
larly at 1.25 pm is scattered in this plane, indicating that 
in these regions the NIR emission from the star forming 
population is lower and the emission traces the evolved 
stellar populations which varies from galaxy to galaxy as 
expected from the different star formation histories. We 
calculate the intensity ratios at 1.25, 3.6, 4.5, 5.6, 8.0 
and 24 pm relative to In a ,int m Table [31 In the right- 
most column, we give the range in possible values for the 
proportionality constant between lu a ,int and intensity at 
other wavelengths within the rms of the fit. 

These results support the conclusion that very red NIR 
colors are due to sources of NIR emission. But the pos- 
sibility that extinction may also be important needs fur- 
ther exploration. In particular, the trend fo r increasingly 
red N IR colors falls along the dust law of Chariot & Fall 
( 2000 ) , which may suggest that the spread in NIR colors 



is at least partly due to dust extinction rather than any 
additional emission mechanism. Using the dust extinc- 
tion image defined by Equation [3j we can measure the 
average dust extinction in the NIR color-selected pixels. 
We further investigate any trends between average inten- 
sities at all observed wavebands and the average dust ex- 
tinction in Figure [5] Again, the data show a correlation 
with the dust extinction. At 5.6, 8.0 and 24 pm larger ex- 
tinction values correlate with higher intensities resulting 
from increased dust emission. These regions, as shown in 
Figure [4j correspond to more intense star forming regions 
and as a result have the most extreme extinction values. 
However, we once again see that the correlation with the 
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Fig. 4. — Surface brightness measurements at Ha, 1.25, 3.6, 4.5, 5.6, 8.0 and 24. /im as a function of the intrinsic Ha intensity (a direct 
measure of the average star formation rate per unit area). Colors are keyed to increasingly red NIR colors (as defined in Table |2l. Linear 
fits with zero offsets are fit to each distribution, whose slope is given in each plot. The rms of the fit is shown by the dotted curves. Only 
points to the right of the dashed vertical line are used for the fit to minimize the contribution of emission due to underlying evolved stellar 
populations. All points are used for the Ha and 24 fim fits as the Ha image is stellar subtracted and the stellar contribution at 24 /im image 
is assumed be negligible. The linear fit to each relation (ie. the intensity ratio relative to lua,int) is tabulated in Table l3j Galaxies with 
anomalous NIR colors (discussed in §|5.2|> whose pixels correspond to compact regions of star formation are shown as purple dots. 



dust extinction is weaker in all wavebands compared with 
the correlation between intensities and intrinsic Ha emis- 
sion shown in Figure [4j In conclusion, Figures [3] to [5| 
indicate that the the NIR color excess is due to emission 
and not dust extinction at wavelengths longer than 1 fim. 

5.2. Galaxies with anomalous NIR color structures 

Not all galaxies show binned pixel distributions sim- 
ilar to the typical ones presented in Figure [2] Ten 
galaxies (15%) show distributions with markedly differ- 
ent structures. To study these, we manually selected 
regions in the NIR color space to identify the spatial 
location of anomalous populations. Their location in 
pixel space and in the spatial domain are shown in pur- 
ple for a 6/10 of these galaxies in Figure [6] (the rest 
can be found in the online edition of Figure ^). For 
these 10 galaxies (Holmberg II, NGC 1097, NGC 1266, 
NGC 1377, NGC 2403, NGC 3031, NGC 4736, NGC 5408, 
NGC 7552, Tololo89) we find that the red colors are re- 
lated to regions of clumpy star formation as indicated 
by compact 24 /am emission. In most of these cases, 
the NIR colors follow a vector corresponding to higher 
I4.5/I3.6 ratios th an the A -0 7 extinction law from [Char- 



iot fc Fall ( 2000 ). For example, Holmberg II, NGC 1266, 
NGC 1377, NUC 5408, NGC 7552 and Tololo89 show a 
similar distribution of NIR colors. The increasingly red 
NIR colors correspond to regions of increasing star for- 
mation (4th panel in Figure [6| and higher dust extinc- 
tion (right panel). These regions correspond to compact 
regions with bright infrared emission. The emission is 
likely related to star formation; in many cases it corre- 
sponds to the brightest regions along a spiral arm or com- 
pact opaque star formation, as is the case for NGC 1377 
which is an opaque starburst galaxy whose dust emission 



is optically thick (Roussel et al. 2006). In this galaxy, 
the emission is related to very young (< lMyr), highly 
obscured starburst activity. These NIR colors are the 
result of either a different extinction law that manifests 
as an optically thick dust screen or perhaps a different 
dust emission law resulting from the high radiation fields 
produced by the starburst activity (in which the ratio 
of PAH emission and continuum emission from hot and 
warm dust is different than the typical spectrum found in 
normal star forming environments found in this study). 
Some a dditional possible mechan isms have recently put 
forth in |Smith fc Hancock ( 2009[ ) to explain the redder 
I4 5/I3 6 ratios seen in nearby dwarf galaxies. 

NGC 3031 and NGC 4736 display NIR color sequences 
that show increasing I4.5/I3.6 ratios for decreasing 
I3.6/I1.25 ratios. Both of these galaxies show sp ectral 
indicators of harboring an AGN ( |Dale et aL]2006| ). 

We have included the points for these 10 galaxies in 
Figures [4j [3] & [5] as purple dots to compare the average 
intensities in these anomalous red colors sequences to the 
star formation rate intensity and dust extinction and find 
that in all cases the values agree with the measurements 
made for the standard NIR color-selected regions. 

6. DISCUSSION 

Spatially resolved pixel-by-pixel analysis of the SINGS 
galaxies yields a number of insights regarding the under- 
lying stellar populations and their associated dust emis- 
sion and extinction. In regions of negligible star forma- 
tion, NIR colors occupy a small region of NIR color space, 
and we have shown that the majority of pixels tend to 
be located in this color space. However, younger pop- 
ulations, with higher Ha and 24 /am emission, tend to 
show colors that become increasingly red in the NIR, 
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Fig. 5. — Surface brightness measurements at Ha (0.656), J (1.25), 3.6, 4.5, 5.6, 8.0 and 24 /im in relation to the average dust extinction 
derived from the ratio between the Ha and 24um intensities as described by Eq.[3] The different colored points are keyed to the regions of 
increasingly larg er N IR colors as defined in §g]and plotted in Figure [2] The purple dots are regions with anomalous NIR color sequences 
as discussed in 515.21 



with both the intensity ratios of I3.6/I1.25 an d I4.5/I3.6 
increasing along a vector well matched b y the A -0 7 ex- 
tinction law from [CharbF^^ll| ([2000|. Motivated by 
this empirical trend in the data, we divided galaxies into 
different regions based on the NIR color of each indi- 
vidual pixel. These regions are designed to encompass 
0.1 dex (or 0.25 mags) above and below the A -0,7 extinc- 
tion curve separated into steps of 0.2 dex (0.5 mags) in 
logI3.6A1.25 color starting at logI 3 . 6 /Ii.25 = -0.6. By 
doing so we isolate regions with different NIR emission 
properties and consequently different stellar populations, 
revealing the relatively linear nature of emission in the 
near- and mid-infrared. 

6.1. The linear scaling of the near /mid- infrared 
spectrum 

Figure [4] shows that the average intensity of NIR color- 
selected regions is strongly correlated with with the aver- 
age intrinsic Ha intensity in that region for the majority 
of bands in this study. The 5.8, 8.0 and 24 jam intensities 
scale linearly with this parameter, with scatter presum- 
ably mostly due the non-negligible contribution from the 
stellar population at 5. 8 and 8.0 /im and al so due to dif- 
ferences in met allicity (Calzetti et al. 12007) and possibly 
PAH fractions ( |Draine"et al.||2007[ ). Rather surprisingly, 
the shorter NIR wavelengths also show a linear correla- 
tion with the star formation rate intensity at sufficiently 
high intensities that the emission is dominated by emis- 
sion from the dust associated with young stellar pop- 
ulations. We find that linear fits to the data with zero 
offsets can fit the data in all wavebands, but limit the fits 
to large star formation intensities (Isfr >10 MJy/sr) in 
the shorter bands and slightly lower star formation in- 
tensities (Isfr >1 MJy/sr) at 5.8 and 8.0 /im. 

The flux ratios between Ha and other wavelengths is 
obtained from the linear fits (and resulting fitting er- 



rors) shown in Figure [4] and tabulated in Table [3j The 
SED points obtained m this way are then normalized 
to the J-band (at 1.25 /im) and plotted as black dots 
in Figure [7J The black open squares represent the col- 
ors (logio I3.eA1.25 = -0.30 ± 0.07 and logio I4.5/I3.6 = 
—0.19 ±0.02) of the majority of pixels in the ensemble of 
galaxies. In Figure [2j this corresponds to the peak NIR 
color in the distribution panel (middle panel) and repre- 
sents the NIR colors for non star-forming environments 
consisting of intermediate and evolved aged stellar popu- 
lations. These NIR colors are consistent with the average 
typical NIR colors fou nd for the integra ted photometry 
of the SINGS sample ( |Dale et aL||2007| ), as well as for 
evolved galaxies in ou r previously studied high-z sample 
(iMentuch et al. 2009). For reference, three SEDs from 



stellar population models using PEGASE.2 are plotted 
for a galaxy with an exponentially declining star forma- 
tion history with an e- folding time of 500 Myr at ages 
of 25 Myr (blue and green lines) and 3 Gyr (red line). 
The green spectrum is for a 25 Myr old galaxy with dust 
extinction in the visible of A v = 0.75, while the blue 
spectrum has no dust extinction applied. These galaxies 
are dominated by the emission from the youngest stellar 
populations and is, in essence, a starburst stellar popula- 
tions modeled according to stellar population synthesis. 
Contribution from non-stellar emission, which is known 
to be signflcant for starburst galaxies at NIR and MIR 
wavelengths is not included. The older spectrum (red) is 
modeled with 0.25 mags of visual extinction. As empha- 
sized earlier (in ^2| , the stellar emission beyond 1 /im is 
not dependent onthe star formation history of the galaxy 
and near-infrared colors for a galaxy that do not have 
emission from PAH and/or dust molecules are indepen- 
dent of a galaxy's age. 

The average spectrum effectively measures the broad 
spectral energy distribution of a young stellar popula- 
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Fig. 6. — IRAC 3.6 /im images for galaxies with anomalous NIR pixel colors are shown for (from top to bottom) NGC 1097, NGC 1266, 
NGC 1377, NGC 2403, NGC 5408 and NGC 7552. The second panel from the left shows the spatial location of the anomalous pixel colors 
(outlined in purple). In general, these colors spatially correspond to bright compact regions of star formation. The left panel of the three 
pixel color panels shows the number of pixels in each color-color bin. The middle panel indicates the median value of the intrinsic Ha 
intensity for all the pixels within that color-color bin. The left panel indicates the median dust extinction for the pixels. 
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TABLE 3 
Linear scaling of Ia with 

^-Ha,int 



X (/im) 


slope med 


an slope 


range 


0.66 


0.22 


0.05 - 


0.55 


1.25 


0.06 


0.02 - 


0.14 


3.60 


0.06 


0.02 - 


0.14 


4.50 


0.05 


0.02 - 


0.10 


5.60 


0.28 


0.11 - 


0.61 


8.00 


0.73 


0.29 - 


1.46 


24.00 


1.07 


0.44 - 


1.72 



tion and its associated dusty star-forming environment 
as the fits were made to minimize the contribution from 
emission from evolved stellar populations (except at Ha 
where all stellar emission is effectively removed). It is 
evident that beyond ~ 2 /im the stellar SED does not 
provide sufficient emission to account for the observa- 
tions. Although SPS models are extremely uncertain 
at these wavelengths, due to poor knowledge of stellar 
evolutionary phases and lack of good spectral templates, 
this conclusion is reinforced by the fact that typical NIR 
colors for stellar populations not contaminated by star- 
forming environments are way lower than those which 
are found in star- forming environments. Emission due 
to dust heated by the young stellar population is the ob- 
vious culprit. We plot for refere nce dust models, s caled 
to the 8.0 /im flux density, from ( |Draine fc Li||2007| ) for 
three PAH fractions of q PA u =0.477o, 2.50% and 4.58% 
(thes e models assume < U > m in = 3.00 and u mSiX = 1000, 
see |Draine fc Li (2007) for details, and are shown for ref- 
erence, not as best tits of the average SED). 

The orange dashed curve is the sum of the PEGASE.2 
star- forming model (blue line) scaled to the 1.25 /im flux 
density and the dust model with #pah = 4.58% scaled to 
the 8.0 /im flux density. The open triangles represent the 
flux densities of this composite model. The model flux 
densities are consistent with the average SED derived 
from the linear fits measured in this analysis. The scaling 
of the near/mid-infrared spectrum with the star forma- 
tion rate suggests this region of the SED for a galaxy can 
be incorporated into galaxy SED models as a spectral 
component that scales linearly with the star formation 
rate. 

6.2. The NIR excess due to circumstellar (disk) 
emission 

A hypothesis for the origin o f the NIR emis s ion we 
have identified was proposed in |Mentuch et al. ([2009 ) , 
namely that the emission is associated with circumstel- 
lar emission of dust around massive young stars. We 
can measure the excess NIR emission at 3.6 /im by sub- 
tracting the contribution of stellar light from the band- 
pass. We use the ratio between emission in the near-IR 
of logI 3 . 6 /Ii.25 = -0-3 (I3.6/I1.25 = 0-5) to measure the 
excess emission in the NIR and plot this as a function 
of the intrinsic Ha. A linear relation exists between the 
NIR excess (stellar subtracted) emission and is shown in 
Figure 8j A linear fit to the data is plotted as a solid 
line. The fit has a slope of 0.0435 (+0.1092/-0.0311). 

We can use the qua ntitative model we developed in 
|Mentuch et al. ( |2009[ ) to see if the excess emission at 



3.6 /im can be explained by the same model. A relation- 
ship between NIR emission from circu mstellar disks from 
the model of |Dullemond et a~ (2001) was applied to the 
integrated light oi a galaxy as: 



V NIR 



(SFR) 



disk 



35OL /M C 



SFR 

M /yr 



yr 



(4) 



this equation requires a SFR in M /yr and a value for 
the timescale of the circumstellar disk emission, t excess . 
The calculatio n assumed an Initial Mass Function from 
Kroupa| fl2001| and as a result requires the same IMF 
assumption when converting the intrinsic Ha intensity 
into a SFR. 



We use the standard conversion from Kennicutt (1998), 
except an additional f actor of 0.55 is applied to con- 
vert from an IMF from |Salpeter (1955) to an IMF from 
Kroupa| ( |200lj ). 



SFR(M /yr) = 0.55 x 7.942 x 10 



-42 



ergs 1 L(H a>int ) 
(5) 



The relation resulting from this model is plotted in Fig- 
ure [8] for two timescales of 10 5 and 10 6 years, which cover 
the range of timescales in which circumstellar disks are 
found locally to survive. In general the model over pre- 
dicts the emission, but is within an order of magnitude 
which is encouraging given the model's simplicity. This 
suggests that the excess NIR emission seen in nearby 
young massive stars can explain the NIR colors seen in 
star forming regions in the SINGS galaxies and may be 
due to circumstellar dust heated to its sublimation tem- 
perature. 

6.3. AGN activity 

Do active galactic nuclei (AGN) cause red NIR col- 
ors? There is a trend for the NIR excess emission re- 
gions to be associated with regions of star formation. For 
many of the normal spiral galaxies in the sample, the re- 
gions of NIR excess spatially correlate with Ha bright re- 
gions, found to be located along the spiral arms. Similar 
knots are seen in the irregular galaxies in the sample, but 
are found dispersed throughout these galaxies. In some 
galaxies (both irregular and spiral types), some of these 
regions are found in the nuclear region of the galaxy. 

Table □ contains a column which identifies whether a 
galaxy may har bor an active gala ctic nucleus using spec- 
tral indicators ( |Dale et aL]|2006[ ). We investigate each 
of the 24 galaxies with an AGN and find that in a third 
of the cases, the nuclear regions of galaxies do corre- 
spond to redder NIR colors. For 5 galaxies (NGC 1266, 
NGC3031, NGC 3198, NGC 4736 and NGC 7552), their 
nuclear regions have very red NIR colors (I3.6/I1.25 > 1 
). The pixels in these regions have both high Hai nt emis- 
sion and high dust extinction. The red NIR colors may 
either be related to the AGN itself (or dust heated by 
the AGN); or could also be from optically thick dust ex- 
tinction in starburst regions as the NIR colors in these 
nuclear regio ns are similar to NGC 1366, an opaque star- 
burst galaxy ( Roussel et al.|[2006 ). 
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Fig. 7. — Using the linear fits from Figure [4] (which scale with the star formation rate), a spectrum of broadband near/mid-infrared 
emission is plotted as solid black points (connected by the black dotted line), normalized at 1.25 jum. This average SED is a measure of the 
young stellar population and its associated dusty star-forming environment as seen in the near/mid-infrared as it is shown to scale with the 
intrinsic Ha emission. The black open squares represent the colors (logio I3.6/I1.25 = —0.30 ± 0.07 and logio I4.5/I3.6 — —0.19 db 0.02) of 
the majority of pixels in a galaxy (taken from the r esults of Figure [2l . S hown for comparison are SEDs generated for stellar emission from 
PEGASE. 2 models and dust emission models from |Draine & Li (2fl 0"7j). The stellar mo dels from PEGASE. 2 are normalized at 1.25 /xm, 
while the dust models, varying in PAH mass fraction (qpAn)i Irom Draine & Li (2007) have been normalized at the IRAC 8 /im band. 
The sum of the young PEGASE model (blue) and qp\n = 4.58% (cyan) dust / PAH model is shown as the orange dashed curve. The open 
triangles represent the flux densities of this model in the 2MASS and Spitzer /IRAC filters and are consistent within errors to the observed 
average SED data points (black solid dots). 



7. CONCLUSION 

The near-infrared, at 1-5 /im, marks a transition in the 
spectral energy distribution of galaxies. At ~l-2 /im, 
emission from longer-lived low-mass stars from evolved 
populations remains bright, tracing the build up of mass 
in a galaxy. However, these wavelengths also trace cur- 
rently star forming populations, which emit brightly in 
both the UV, optical and near-infrared. Near-infrared 
colors at 1-5 /im, if dominated by stellar emission, are 
predicted to be independent of the age of a stellar popula- 
tion. However, in star-forming environments, the contri- 
bution from Polycyclic Aromatic Hydrocarbons and dust 
grain emission becomes important beyond 2 /im, making 
the near- and mid-infrared a complicated but rich region 
of a galaxy's spectral energy distribution. 

In this paper, we have analyzed the pixel- by-pixel near- 
infrared colors for a sample of 68 galaxies from the Spitzer 



Infrar ed Nearby Galaxy survey (SINGS; Kennicutt et al. 
2003) and Large Galaxy Atlas (LGA; |Jarrett et al |2UU3p . 
Specifically we use images at 1.25, 3.6, 4.5, 5.6, 8.U and 
24 /im, as well as continuum subtracted Ha images for 
56 galaxies. Using the Ha and 24 jam images, we were 
able to create images representing the spatially resolved 
dust extinction (at A(Ha)) and the intrinsic Ha inten- 
sity. The pixels in each galaxy were categorized into 
different regions based on a near-infrared color selection. 



The majority of pixels for each galaxy in the sample are 
found to have constant NIR colors, as is predicted for 
main-sequence stellar emission, while the reddest NIR 
pixels correspond to regions of star formation, which usu- 
ally trace the galaxy's spiral arms structure. Our anal- 
ysis suggests that very red NIR colors are produced by 
emission, not dust extinction. 

For each NIR color cut in our galaxy sample, we mea- 
sured the average intensities at 1.25, 3.6, 4.5, 5.6, 8.0 
and 24 /im for each galaxy. We compared these to the 
average dust extinction and star formation rate in each 
region and found a linear correlation between the emis- 
sion at all bands from ~l-8/im with the average Hai nt 
intensity (and consequently the star formation rate). The 
strong correlation between near-infrared and PAH emis- 
sion from 1-8 /im and star formation suggests the average 
spectrum derived from this dataset is representative of a 
young stellar population and its associated dusty star- 
forming environment and can be incorporated to galaxy 
SED models very simply as a component directly scaled 
to the star formation rate of a stellar population. 



This analysis relied entirely on archival data. We 
would like t o thank the SINGS collaboration (Kennicutt 
et al. 2003[ ) and the Large Galaxy Atlas team jjarrett 




Ha,int 

Fig. 8. — Excess NIR emission is plotted by subtracting the con- 
tribution of stellar emission at 3.6 /im using emission at 1.25 /im. 
The excess emission is shown to correlate with the star formation 
rate with a median slope of 0.0435(+0.1092/-0.0311). The dashed 
and dotted lines is the equatio n for NIR exces s emission due to 
circumstellar disk emission from [Mentuch eT al. (2009), as applied 
m Eq. Eland [SI for disk lifetimes of 10 b and 10 13 years respectively. 
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